Finite width and local field corrections to spin Coulomb drag in a 
quasi-two-dimensional electron gas 
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We study the spin Coulomb drag in a quasi-two-dimensional electron gas of finite transverse width, 
including local field corrections beyond the random phase approximation (RPA). We find that the 
finite transverse width of the electron gas causes a significant reduction of the spin Coulomb drag. 
This reduction, however, is largely compensated by the enhancement coming from the inclusion of 
many-body local field effects beyond the RPA, thereby restoring good agreement with the experi- 
mental observations by C. P. Weber et al, Nature, 437, 1330 (2005). 
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Lately a number of spin-based devices have been de- 
signed and studied with the goal of combining memory 
and logical functions in a single device [E 0i S 0, 01 ■ ^ 
crucial requirement for most spin devices is the existence 
of a robust population of spin-polarized carriers Q as 
well as the feasibility of moving the spins in the device 
by means of spin currents. In this context, understanding 
the various scattering processes which control the relax- 
ation of the spin and the spin current has emerged as 
an important theoretical and practical problem. In this 
Letter we focus on the spin current relaxation caused 
by electron-electron scattering in quasi-two-dimensional 
quantum wells of finite width. 

It is now well established [1, S, 0| that the Coulomb 
interaction induces momentum transfer between the car- 
rier populations of "up" and "down" spin. This results 
in a spin Coulomb drag (SCD) effect within a single two- 
dimensional electron gas (2DEG) layer, in close analogy 
to the conventional charge Coulomb drag, which occurs 
between spatially separated layers. Because the inherent 
friction between the two spin components leads to a decay 
of the spin current even in the absence of impurities, the 
SCD has recently become a subject of exp erimental and 



theoretical investigations [7|, [10|, |11 , 13, ij- In particular, 
the SCD is known to reduce the spin diffusion constant 
relative to the conventional density diffusion constant, 
and thus to prolong the time during which a spin packet 
can be effectively manipulated. Indeed a significant re- 
duction in the spin diffusion constant was measured by 
Weber et al. and found to be in quantitative agree- 
ment with the SCD reduction calculated for a strictly 
2DEG in the random phase approximation (RPA) [loj . 

The agreement between the RPA theory for an electron 
gas of zero width and the experiment is encouraging, but 
also somewhat puzzling. Quantum wells in which the 
experiments have been performed are not strictly two- 
dimensional but have finite transverse widths of the or- 



der of 10 nm. In contrast to the conventional charge 
Coulomb drag, where the inter-layer spacing causes an 
exponential suppression of large angle scattering events, 
the main contribution to the SCD comes from events with 
a momentum transfer of the order of the Fermi momen- 
tum. At such values of the momentum transfer, the form 
factor that takes into account the width of the 2DEG is 
significantly smaller than 1, and should cause a signifi- 
cant reduction of the SCD even at relatively high tem- 
perature and densities. 

A second problem is that the RPA on which most cal- 
culations of SCD [l3| so far have been based, ne- 
glects some important physical effects. Namely, in RPA 
one ignores the fact that the presence of an electron at a 
point in space reduces the probability of other electrons 
coming nearby, via the so-called exchange-correlation 
(xc) hole 15|. The repercussions of this effect upon the 



effective electron-electron interaction are quite subtle be- 
cause the xc hole, while keeping the electrons apart, also 
reduces the collective screening of the interaction. Re- 
cently we have studied the effect of the xc hole in charge 
Coulomb drag and found that its introduction results in a 
significant enhancement of the charge transresistivity and 
qualitative changes in its temperature dependence 
The sheer size of these effects suggests that a quantitative 
description of the SCD cannot be achieved by a theory 
that does not include xc effects. 

In this Letter we report the results of our calculations 
of the SCD taking into account (i) the finite transverse 
width of the quantum well and (ii) xc effects beyond the 
RPA [ii], Following Ref. [ij, we have used many- 

body local field factors [l5| to take into account the modi- 



fications to the effective electron-electron interaction [16 1 
due to the xc hole. The main outcome of our study is 
that the enhancement of spin drag caused by the many- 
body local field factors largely compensates the reduction 
of the effect coming from the finite well width, thereby 
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restoring good agreement with the experimental obser- 
vations of Ref. 0- 

SCD manifests itself when a current of active spin- 
down electrons, J|, through the Coulomb interaction, 
transfers momentum to the passive spin-up electrons with 
zero current, J| = 0. This process induces a gradient of 
electrochemical potential, which results in a "spin elec- 
tric field" The direct measure of SCD is the spin 
drag resistivity, defined as 



PliiT) ^E-il'h when J| = 



(1) 



We calculate the spin drag resistivity in the ballistic 
regime as appropriate for the clean samples used in the 
experiments '7'|. The spin transresistivity is negative. Its 
absolute value at temperature T is given by [a| 







Imlloi (q, cj)Imnox (g, uj) 
sinh^(fiw/2T) 



(2) 



where huj and hlf are, respectively, the transferred en- 
ergy and in-plane momentum from the spin a ='\ to the 
spin a =i, Wi^i{q,Lo) is the dynamically screened effec- 
tive interaction between electrons with spins t and |, 
IlQa{q,uj) is the spin-cr component of the non-interacting 
finite temperature polarization function 15|, n| and rii 



are the densities of up and down spins respectively, and 
A is the area of the sample. Further we will assume 
that the electron gas has no net spin polarization, i.e. 

To calculate the effective interaction W||(a;) we have 
used an approximation scheme developed a few years 
ago by Vignale and Singwi (VS) following the_ 
oneering treatment of Kukkonen and Overhauser 
Within this scheme the effective electron-electron inter- 
action, W{q,uj) (a 2 X 2 symmetric matrix with compo- 
nents W||, W^i, and W^i) is approximated as 

W{q, Lu) = viq) ~ Viq, Lu) ■ fl{q, to) ■ V{q, Lo) (3) 

where the full polarization function, n(g,a;), is defined 
in terms of the non-interacting polarization function, 
no(q,w), and the unscreened effective Coulomb interac- 
tions, V {q, Lo) (defined below) has the form 



%,c^)=no(g,a;)- I + V{q,Lo) ■Mq,uj) 



(4) 



Here all quantities are 2x2 matrices and the dot denotes 
a matrix product. Selecting the f i component we get 

Wuil, = V^Ti -I- v{q)G^i (g, uj)F{qd) , 

(5) 

where the unscreened effective Coulomb interaction be- 
tween spins a and a' is 

V„„, {q, lu) = v{q) (1 - (q, u;)) F{qd) , (6) 



where d is the width of the well. Observe how the un- 
screened effective interaction Vaa'iq,^) differs from the 
bare Coulomb interaction v{q) ~ 27re^/Kog in two ways. 
First, the spin-resolved local field factors Go-o-' reduce 
the bare interaction by a factor of 1 — G^^i {q,uj): this 
is the effect of the xc hole, mentioned in the introduc- 
tion. Second, the form factor F{qd) takes into account 
the form of the wave function of the lowest transverse 
subband, and is defined as the square of the Fourier trans- 
form of the corresponding density profile p{z) (for p{z) = 
{2/d) sin(7rz/d)2 we have F{qd) = I - (1/3 - 5/47r2)gd 
for qd ^ and F{qd) 3/(47r^qd) for qd ^ oo). Taking 
d = 12 nm and q c:^ kp we have ~ 2 for n = 4.3 ■ 10^^ 
cm~2 so that \F{qd)\'^ ~ 0.5. 

In Eq. ([5]) e{q,u!) is an effective dielectric function, 
which can be represented in the following manner 

e{q,u) = {l + V+iq,Lo)Uoiq,Lo))il + V-iq,Lo)Uo{q,Lu)) , 

where V±{q,uj) = {V^^{q,uj) ±V-^i{q,uj)) /2. Introduc- 
ing the corresponding notation for the local field factors, 
G± = (G|| ± G||)/2, we can write 

V+{q,L0) = v{q)il~G+iq,Lu))F{qd) , (7) 
V^{q,i0) = -v{q)G-{q,Lu)Fiqd) (8) 

G+ and G_ are known as the "charge-channel" and the 
"spin-channel" local field factors respectively. 

In our actual calculations, described below, the lo- 
cal field factors are derived from diffusion Monte Carlo 
simulations 1^ 1^) in the low-energy finite-momentum 
transfer regime, and from a combination of analytical and 
empirical methods [20l.[2]| in the small- momentum, high- 
energy transfer regime. These two regimes correspond to 
two distinct mechanism of electron-electron interaction: 
in the first, the interaction is mediated by electron-hole 
pairs, and in the second by plasmons ^23\. 

Notice that both the dynamic and static local field fac- 
tors 1^, 21 1, which we have used in our calculations, are 
evaluated at T = for an ideal 2DEG of zero width. At 
present there are no calculations of the local field factors 
that take into account the finite width of the quantum 
well and/or the temperature dependence. We expect, 
however, that the temperature dependence of the local 
field factors will play a relatively minor role in compar- 
ison to the temperature dependence that we explicitly 
include in the Fermi factors and in the non-interacting 
polarization functions. Similarly, the effect of the finite 
width of the well should be largely taken care of by our 
use of the form factors F{qd) in Eqs. ^ and ([6|). 

Recently Weber et al. [7] have demonstrated experi- 
mentally that the spin diffusion constant of electrons in 
an n-type GaAs quantum well is strongly reduced by the 
Coulomb interaction over a broad range of temperatures 
and electron densities. Their measurements were in good 
agreement with the theoretical calculations of D'Amico 
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FIG. 1: Spin drag resistivity p^i as a function of temperature, 
calculated within RPA. The solid, dashed, and dotted curves 
correspond, respectively, to pn for d = 0, 12, and 20 nm. The 
black squares represent p^i, deduced from the experimental 
data reported in Fig. 3 of Ref. for the sample with n = 
4.3 ■ 10^^ cm~^ and d = 12 nm. 
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FIG. 2: The scaled spin drag resistivity p^i/pD vs T for 
= 1. The solid and dash-dotted curves are calculated for a 
2DEG of the width d = Q and 10 nm within the VS scheme. 
The dashed curve corresponds to P]i/pD of an ideal 2DEG, 
calculated within the RPA. The ordinary Drude resistivity, 
pD, is calculated for the mobility p — 300 V cm^^ s^^. 



and Vignale [13], which, however, neglected both the fi- 
nite width of the quantum well and the local field factor. 
Here we present results of calculations which include both 
these effects. 

First we show that the effect of the finite width of 
the quantum well is very important. In Fig. [1] we plot 
the absolute value of the spin drag resistivity, calculated 
within RPA, vs temperature for three different values of 
the width, d = 0, 12, and 20 nm. The solid curve repro- 
duces the RPA calculations of the spin drag resistivity 
for an ideal 2DEG of zero width 10]. The black squares 
show the spin drag resistivity, as deduced from the exper- 
imentally measured spin diffusion constant Ds according 
to the formula 0, HI] 

Here the charge diffusion constant Dc is obtained from 
the electric resistivity p via the Einstein relation, and 
Hs/Hq is the many-body enhancement of the spin sus- 
ceptibility of the electron gas, which, for the time being, 
following Ref. 0, we approximate as 1. Fig. [T] shows that 
an increase in the width of the quantum well is accom- 
panied by a strong reduction of the spin drag resistivity: 
aiT ^Tf {Tf ~ 178 K) p|| for d = 12 nm is a factor 
of two smaller than at d = 0. The physical reason for 
this effect is that the SCD is dominated by large angle 
scattering events, where the momentum transfer q is of 
the order of the diameter of the Fermi sphere. At these 
values of q the form factor F{qd) is substantially smaller 
than 1 - the more so the higher the density. By contrast, 
in the ordinary charge Coulomb drag large angle scat- 
tering events are exponentially suppressed for g > 1/A 
where A is the inter-layer separation. Then, the finite 
width of the quantum well produces corrections of the 
order of \F{d/A)\'^, which is rather close to 1. 



Next we include the effect of exchange and correlation 
on the spin drag resistivity. In Fig. [5] we compare the 
spin drag resistivity, calculated within the VS scheme, 
with the corresponding result obtained within the RPA. 
The scaled spin drag resistivity is shown as a function of 
T for the dimensionless parameter = 1 corresponding 
to ri = 3.5 • 10^^ cm^^ and to the Fermi temperature 
Tf — 147 K. Notice that the inclusion of the local field 
factors shifts the peak of the spin drag resistivity towards 
lower temperatures compared to the RPA. Accordingly, 
the spin drag resistivity curves show quite different slopes 
at high temperatures depending on whether local field 
factors are included or not. As a result, the reduction 
in spin drag coming from the finite width of the well 
is largely compensated by the many-body enhancement 
of the spin drag resistivity. The situation is similar for 
= 2. We conclude that the combined effect of the 
finite width and the xc hole on the effective interaction 
restores the agreement between theory and experiment. 

In Fig. [3] we verify this conclusion for the strict ex- 
perimental situation of Ref. 0- Here the dotted and 
dashed curves are calculated within RPA, respectively, 
for 0? = and 12 nm and show the strong reduction of 
spin drag due to the form factor of the quantum well. 
The spin drag resistivity is deduced from the experi- 
mental data for the spin diffusion via Eq. ([9]). How- 
ever, we have now taken into account the many-body 
enhancement of the spin susceptibility with respect to 
the non-interacting spin susceptibility. In a paramag- 
netic electron liquid the Hartree terms cancel and the 
enhancement is determined by the "spin-channel" many- 
body local field factor G-{q, 0) according to the formula 
Ho/Us = 1-l-limg^o V-{q, 0)Ho(g, 0) [H]. Then, it is seen 
that our theoretical curve for the spin drag resistivity, in- 
cluding both the many-body effects beyond RPA and the 
finite width of the quantum well, agrees very well with 
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FIG. 3: Spin drag resistivity as a function of temperature 
for the electron density n = 4.3 ■ lO^'^ cm~^. The dashed 
and dotted curves correspond to pn, calculated within the 
RPA for d = and 12 nm, respectively. The solid curve 
correspond to p-^y, calculated beyond the RPA within the VS 
scheme for d = 12 nm. The symbols represent pn, deduced 
from the experimental data for d = 12 nm of Ref. [Tj. Notice 
that the data points for the spin drag resistivity are lower 
than in Fig. 1, because we have taken into account the many- 
body enhancement of the spin susceptibility in converting the 
experimental values of Ds to according to Eq. ([9]). 



the experimental data of Ref. while the RPA, in the 
presence of the form factor, does not. 
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Our calculations do not include the spin-orbit interac- 
tion in the one-electron states because the measurement 
of Ref. 0] were done on a spin grating with wave vector 
along the [11] direction of GaAs, which is the direction 
along which spin-orbit effects are essentially non existent. 
We might add that even in different directions the influ- 
ence of the spin-orbit coupling on the spin diffusion con- 
stant would show up as a significant effect only at wave 
vectors much smaller than the one studied in Ref. 01 • Fur- 
thermore, according to a recent calculation in Ref. p^ . 
the spin Coulomb drag is enhanced by Rashba spin-orbit 
coupling, 7, by a small factor 1-1-7^, which for the typical 
value of 7 ~ 0.1, is a 1% correction. 

We have neglected the phonon-mediated SCD. In a single 
electron layer the phonon-mediated interaction between 
electrons cannot compete with the Coulomb interaction, 
due to the smallness of the electron-phonon coupling con- 
stant. The phonon-mediated interaction is effective only 
in bilayers with large barrier thickness, where the inter- 
layer Coulomb interaction becomes very small. 
In the static regime the local field factors are deduced 
from the computed static response of 2DEG to periodic 
external potentials of wave vector q: they vary linearly 
with q both at small q and at large q but with different 
slopes as discussed in Section 5.6 of Ref. [l^. The local 
field factor in the dynamical regime is obtained from a 
perturbative mode-coupling theory, described in [20I. [2l|. 



